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The properties of the momentum space quark propagator in Landau gauge are examined for the overlap quark 
action in quenched lattice QCD. Numerical calculations are done on three lattices with different lattice spacings 
and similar physical volumes to explore the approach of the quark propagator towards the continuum limit. We 
have calculated the nonperturbative momentum-dependent wavefunction renormalization function Z{p'^) and the 
nonperturbative mass function M{p^) for a variety of bare quark masses and extrapolate to the chiral limit. We 
find the behavior of Z{p^) and M(p^) are in good agreement for the two finer lattices in the chiral limit. The 
quark condensate is also calculated. 



1. INTRODUCTION 

There have been several studies of the momen- 
tum space quark propagator |1I2I3I4|5] using dif- 
ferent gauge fixing and fermion actions. Here we 
focus on Landau gauge fixing and the overlap- 
fermion action and extend previous work '5' to 
three lattices with different lattice spacings a and 
very similar physical volumes. This allows us to 
probe the scaling behavior and the continuum 
limit of the quark propagator in Landau gauge. 

2. QUARK PROPAGATOR ON THE 
LATTICE 

In the continuum the rcnormalized Euclidean- 
space quark propagator must have the form 



1 



(1) 



where Z{C^;p^) is the wavefunction renormaliza- 
tion function, M{p'^) is the nonperturbative mass 
function, and C is the renormalization point. 

On the lattice the bare quark propagator can 
be written as 



J2c,.{phA + B{p) . 



(2) 
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We use periodic boundary conditions in the spa- 
tial directions and anti-periodic in the time di- 
rection. The discrete momentum values for a 
lattice of size iV^^ x Nt , with rii = 1 , . . , TV, and 
nt = 1, .., iVt, are 



Pi = 



2tv 
Ma 



2 2 



We can perform a spinor and color trace to iden- 
tify 



CM = ^tr[7M5^ 



= S(p) = ^tr[5'^-^(p)]. 



The general approach to tree-level correction|2| 
utilizes the fact that QCD is asymptotically free 
and so it is the difference of bare quantities from 
their tree-level form on the lattice that contains 
the best estimate of the nonperturbative infor- 
mation. For the overlap fermion, the tree-level 
correction is nothing but to identify appropriate 
kinematic lattice momentum q. We can identify 
the appropriate kinematic lattice momentum q di- 
rectly from the definition of the tree-level quark 
propagator numerically. 



Cl"\p) 



(C(0)(p))2 + (S(0)(p))2 



(3) 



We can also obtain the kinematic lattice momen- 
tum q analytically 0]. 
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Having identified the lattice momentum q, we 
can now define tlie bare lattice propagator as 



^barc 



iiA{p) + B{p) 
Z,{C;a)S{C;p).. 



ii + M (p) 



(4) 



where S{(^;p) is the lattice version of the renor- 
malized propagator in Eq. and Z2{C',ci) is 
the quark wave-function renormalization con- 
stant chosen so as to ensure Z{C; C^) = 1. 

The overlap fermion formalism 6 7 realizes an 
exact chiral symmetry on the lattice and is au- 
tomatically 0{a) improved. The massive overlap 
operator can be written as as |H| 



(5) 



where p is the Wilson mass with negative sign, 
and the quark mass parameter 77 = mo/2p. Writ- 
ten according to bare quark mass mg, we have 



D(mo) 



1 

2~p 



P+ — + (P- — )75e(-ff) 



(6) 



and the overlap quark propagator is given by 
S^'^^^mo) ^ D-\r^) , (7) 
where 



2p 



1 



[D-\v)-l] 



1-77 

3. NUMERICAL RESULTS 



(8) 
(9) 



Here we work on three lattices with different 
lattice spacing a and very similar physical vol- 
umes using a tadpole-improved plaquette plus 
rectangle gauge action. For each lattice size, 50 
configurations are used. Lattice parameters are 
summarized in Tabled 

In the calculations, k = 0.19163 was used for 
all three lattices, which gives pa= (8 — 1/k)/2 = 
1.391. We calculate for 10 quark masses on each 
lattice by using a shifted Conjugate Gradient 
solver. The 14th order Zolotarev rational ap- 
proximation is used to evaluate the matrix sign 
function e[Hw). The ten bare quark masses we 



Table 1 

Lattice parameters. 



Action 


Volume 




a (fm) 


Wo 


Improved 


16^ X 32 


4.80 


0.093 


0.89650 


Improved 


12^ X 24 


4.60 


0.125 


0.88888 


Improved 


8^ X 16 


4.286 


0.194 


0.87209 



use in our calculation are toq = 2/3?7 — 106, 124, 
142, 177, 212, 266, 354, 442, 531, and 620 MeV 
respectively. 

The detailed results will be presented else- 
where. Here we focus on the comparison of the 
results on these three lattices. All data has been 
cylinder cut and extrapolated to the chiral 
limit using a simple linear extrapolation. The 
mass function M{p) in the chiral limit for the 
three lattices is plotted in Fig. ^ and the renor- 
malization function Z{p) of the three lattices is 
plotted in Fig.|21 We can see that when the mass 
function M{p) is plotted against the discrete lat- 
tice momentum p the results of the three lattices 
arc in good agreement, while for the renormaliza- 
tion function Z{p), good agreement is reached on 
the three lattices if it is plotted against the kine- 
matical lattice momentum q. The overall agree- 
ment between the two finer lattices is good. 

In the chiral and continuum limits, the asymp- 
totic quark mass function has the form 



[ln(^VA2 



QCD 
QCD 



pi 



(10) 



where the anomalous dimension of the quark mass 
is (Im = 12/(33 — 2Nf). Using the momentum p, 
in the fitting range ap c(1.3, 2.5), on the 16^ x 32 
lattice, the resulting value for the quark conden- 
sate is 



(tPiP) = -(288±24MeV)^ 



(11) 



This is in excellent agreement with the value 
(^■0) = -(268±27MeV)3 extracted from the Asq- 
tad action using the same method jH]]. 

4. SUMMARY 

In this report, we have considered tadpole- 
improved quenched lattice configurations, and the 
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Figure 1. Comparison of tlie mass function M{p) of 
tliree lattices in tlie cliiral limit. The upper graph 
is plotted against the discrete lattice momentum p 
and the lower graph is plotted against the kinematical 
lattice momentum q. 

overlap fermion operator with the Wilson fermion 
kernel. The momentum space quark propagator 
has been calculated in Landau gauge on three lat- 
tices with different lattice spacing a and similar 
physical volumes to explore the scaling property. 
The continuum limit for Z{j)) is most rapidly ap- 
proached when it is plotted against the kinemati- 
cal lattice momentum g, whereas the quark mass 
function, M{p), should be plotted against the dis- 
crete lattice momentum p. The good agreement 
between the two finer lattices suggests that we 
are close to the continuum limit. 
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Figure 2. Comparison the renormalization function 
Z{p) — Z{(^;p) for renormalization point =5.31 GeV 
for the g-scale and =3.15 GeV for the p-scale for the 
three lattices in the chiral limit. The upper graph 
is plotted against the discrete lattice momentum p 
and the lower graph is plotted against the kinematical 
lattice momentum q. 
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